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Nucleon form factors and O(a) Improvement* 
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Nucleon form factors have been extensively studied both experimentally and theoretically for many 
years. We report here on new results of a high statistics quenched lattice QCD calculation of vector and 
axial- vector nucleon form factors at low momentum transfer within the Symanzik improvement programme. 
The simulations are performed at three k and three j3 values allowing first an extrapolation to the chiral 
limit and then an extrapolation in the lattice spacing to the continuum limit. The computations are all 
fully non-perturbative. A comparison with experimental results is made. 



1. INTRODUCTION 

For many years experiments have been per- 
formed with electron-nucleon scattering to ob- 
tain information about the structure of the nu- 
cleon. Form factors are defined from the gen- 
eral decomposition of the proton, p (or neu- 
tron, n) matrix element]^] (q = p — p'): 



Experiments lead to phenomcnological dipole 
fits: 



Gl(-q 2 ) 



G n e (-q 2 



G p m (-q 2 



G n m (-q 2 ) 



= i/(i + H 2 A4)) 2 , 
~ o, 



-1.91. 



u(p, s) [ ltl Fl + a^Fi] utf, s>). 

We have Fi(0) = 1 as V is a conserved current, 
while F 2 (0) = fi — 1 measures the anomalous 
magnetic moment (in magnetons) . Usually we 
define the Sachs form factors: 



Ge(-q 2 ) 



Fi(-q 2 ) + 



F2(-q 2 ), 



(2m) 2 

G m (-q 2 ) = M~q 2 ) + F 2 (-q 2 ). 

*Talk given by R. Horsley at Lat98, Boulder, U.S.A. 
x We have already re- written everything in euclidean 



space, so that eg p = (iE p ,p) and —q 



JM)2 > a 



with my ~ 0.82 GeV, ~ 2.79, fi n < 
Neutrino-neutron scattering, nv^ 
gives from the charged weak current the axial 
form factor g A (—q 2 ). In addition g A (0) is also 
accurately obtained from /?-decay, n — > pe~V. 
Upon using current algebra this form factor 
can be related to the matrix element: 

The phenomcnological fits are: 
g A (q 2 )=g A (Q)/{l+{-q 2 /m\)) 2 , 
with g A (0) = 1.26, m A ~ 1.00 GeV. 
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2. THE LATTICE METHOD 



Quenched configurations have been gener- 
ated at /3 = 6.0 (0(500), 16 3 x 32 lattice) 
= 6.2 (0(300), 24 3 x 48 lattice) and (3 = 6.4 
(0(100), 32 3 x 48 lattice), @]. By forming the 
ratio of three-to-two point functions, |2|: 



(N(t;p)N(0:,p)) 



{N(T:p)N(0;p))(N(t;p)N(0;p))(N(t-T:p')N(0;p')) 
(N(r:p')'N(Q;p')){N(t;p')7l(Q;p')){N(t-T-p)W(0;p)) 



1st Order Pert. Theory 



Tl Pert. Theory 



■ <N|V 4 |N> 

• a result 



cx (N a {p)\0{q)\Np{^)), 

the appropriate matrix elements can be found. 
(Only the quark line connected part of the 
3-point function is considered.) For each /? 
we chose three k values and a variety of 3- 
momenta: p = 2n/N s { (0,0,0), (1,0,0) }, 
q = 2n/N s { (0, 0, 0), (0, 1, 0), (0, 2, 0), (1, 0, 0), 
(2,0,0), (1,1,0), (1,1,1), (0,0,1) } together 
with the nucleon either unpolarised or po- 
larised in the y direction. (Some combinations 
were too noisy to be used though.) After sort- 
ing the matrix elements into q 2 classes (de- 
fined by q 2 in the chiral limit), 4-parameter 
fits are made assuming that the form factors 
are linear in the bare quark mass am q . 0(a) 
improved Symanzik operators are used: 

= Z v (l + b v am q ) x 

A* = Z A (1 + b A am q ) x 

[V'7 M 75'0 + IcAad^^ip)] , 

where Zy, Za, by, cy, ca (and c sw ) have 
been non-perturbatively calculated by the Al- 
pha collaboration, ||. All matrix elements 
thus are correct to 0(a 2 ). We can check Zy 
as V M is a conserved current (ie -Fi(O) = 1). 
In Fig. |l| we show a comparison of the two 
determinations of Zy ■ Very good agreement 
is seen. This is not the case when Wilson 
fermions are used (see ref. g). Finally we 
note that although we have included the im- 
provement terms in our operators, numerically 
they seem to have little influence on the value 
of the matrix element. 



Figure 1. Zy for improved fermions. Shown is 
the lowest order perturbation result together with a 
tadpole-improved version (as given in M). The non- 
perturbative determinations are shown as open circles, 
|3|, and filled squares, this work. 
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Figure 2. The proton form-factors Ge(-q 2 ) and 
Gm(—Q 2 ) against — q 2 showing experimental results 
(open circles, taken from ref. [H) and lattice results 
(filled circles, /3 = 6.2 only). The string tension is 
used to fix the scale as in |4 1. All fits are dipole fits. 



3. RESULTS 

In Fig. | we show G p e {-q 2 ) and G p m {-q 2 ) 
for (3 = 6.2 together with experimental results 
(also plotting the other [3 values tends to clut- 
ter the picture). Making dipole fits gives Fig. || 
for the continuum extrapolation. There seems 
to be little inclination for my to approach the 
experimental result. (A roughly similar result 
is obtained from G^, although due to larger 
error bars the results are more compatible.) 

For the axial current we find the results in 
Figs. ||, [| The form factor fall-off is again too 
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Figure 3. m v from f3 = 6.0, 6.2, 6.4 against a 2 . The 
phenomcnological value is also given at a 2 = 0. 




Figure 4. gA(—g 2 ) against -q 2 , notation as in Fig.]^ 

soft as via is too large. However the impor- 
tant 3,4(0) is faring better, see Fig. ||. 



4. CONCLUSIONS 

We have performed simulations at three (3 
values so that an attempt can be made to 
take the continuum extrapolation, a — > 0. 
While the lattice dipole masses seem to be too 
large, <7a(0) is in reasonable agreement with 
the experimental result. The mass discrep- 
ancies may be due to a quenching effect, al- 
though only similar simulations using dynam- 
ical fermions will be able to answer this. 
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Figure 5. The continuum extrapolation of m,A- 
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Fi glire 6. The continuum extrapolation of <?a(0). 
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